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We describe a general framework to study covariant symmetric broadcasting maps for mixed 
qubit states. We explicitly derive the optimal N — > M superbroadcasting maps, achieving optimal 
purification of the single-site output copy, in both the universal and the phase covariant cases. We 
also study the bipartite entanglement properties of the superbroadcast states. 



I. INTRODUCTION 



"Information" is by its nature broadcastable. What happens when information is quantum and we need to distribute 
it among many users? Indeed, this may be useful in all situations where quantum information is required in sharable 
form, e. g. in distributed quantum computation for quantum shared secrecy Q, and, generally, in quantum game- 
theoretical contexts 3]. However, contrarily to the case of classical information, which can be distributed at will, 
broadcasting quantum information can be done only in a limited fashion. Indeed, for pure states ideal broadcasting 
is equivalent to the so-called "quantum cloning" , which is impossible due to the well-known "no-cloning" theorem 
[HOI El (see also 0,HEJE3)- The situation is more involved when the input states are mixed, since broadcasting can 
be achieved with an output joint state which is indistinguishable from the tensor product of local mixed states from 
the point of view of individual receivers. Therefore, the no cloning theorem cannot logically exclude the possibility of 
ideal broadcasting for sufficiently mixed states. 

In Ref. 1 1 1 1] it was proved that perfect broadcasting is impossible from a single input copy to two output copies 
for an input set of non mutually commuting density operators. This result was then considered (see Refs. 1 1 1 1 ] and 
|12| ) as an evidence of the general impossibility of broadcasting mixed states drawn from a non-commuting set in 
a more general scenario, where N ^ 1 equally prepared input copies are broadcast to M > N users. However, for 
sufficiently many input copies N and sufficiently mixed input states the no-broadcasting theorem does not generally 
hold |l3j ] , and for input mixed states drawn from a noncommuting set it is possible to generate M > N output local 
mixed states which are identical to the input ones, by a joint correlated state. Actually, as proved in Ref. 13], it is 
even possible to partially purify the local state in the broadcasting process, for sufficiently mixed input states. Such 
a process of simultaneous purification and broadcasting was named superbroadcasting. For qubits, the fully covariant 
superbroadcasting channel that maximizes the output purity (i. e. the length of the output Bloch vectors of local 
states) when applied to input pure states coincides with the optimal cloning map |14| . 

The possibility of superbroadcasting does not increase the available information about the original input state, due 
to unavoidable detrimental correlations among the local broadcast copies, which do not allow to exploit their statistics 
(a similar phenomenon was already noticed in Ref. |l5j|). Essentially, superbroadcasting transfers noise from local 
states to correlations. From the point of view of single users, however, the protocol is a purification in all respects, 
and this opens new interesting perspectives in the ability of distributing quantum information in a noisy environment, 
and deserves to be analyzed in de pth . For qubits, it has been shown that for universal superbroadcasting is possible 
with at least N = 4 input copies |l3j . Is this the absolute minimum number for superbroadcasting, or does it hold 
only for this particular set of input states? In this paper we will show that, indeed, for equatorial mixed states of 
qubits the minimum number is N = 3. However, for smaller non-commuting sets of qubit states the possibility of 
superbroadcasting with only N = 2 input co pies is still an open problem (for larger dimension d > 2 it is possible to 
superbroadcast also for N = 2, see e.g. Ref. flq). 

We want to point put that clearly there are limitations to superbroadcasting. The input state must be indeed 
sufficiently mixed, since pure states cannot be broadcast by the no cloning theorem. However, states with a pretty 
high purity can sill be superbroadcast, e. g. for universally covariant superbroadcasting |13| from N = 4 to M = 5 it 
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is possible to superbroadcast states with a Bloch vector length up to 0.787 (0.935 for the phase covariant case). One 
can achieve superbroadcasting with even higher input purity for increasing N, approaching unit input Bloch vector 
length in the limit of infinitely many input copies. There are also some limitations in the absolute number of output 
copies for which one can achieve superbroadcasting. E. g. in the universal case, for N = 4 input copies one can 
superbroadcast up to M — 7 output copies, for N = 5 up to M — 22, and for N > 5 up to infinitely many. The 
output purity is clearly decreasing versus M. In this paper we will further analyze all limitations to superbroadcasting 
also for the case of equatorial input qubits. 

Regarding the possibility of achieving superbroadcasting experimentally, the first route to explore is to use the 
same techniques as for purification , since the superbroadcasting map for M > N generalizes the purification map 
using the same protocol |l8j . This transformation involves a measurement of the total angular momentum of the 
qubits, then an optimal Werner cloning [l4| in the universal case (or an optimal phase covariant cloning [24| in the 
phase covariant case) . Another possibility is to use the methods of Ref. in order to classify all possible unitary 
realizations, and then seek for experimentally achievable ones using current technology. 

In this paper we present phase covariant superbroadcasting, and also give a complete derivation of the universal 
superbroadcasting map, presented in Ref. [Tflj . The two maps are derived in a unified theoretical framework. The 
paper is organized as follows. In Section [H] we introduce some preliminary notions regarding symmetric covariant 
maps. In Sects. IIIII and llVI we give a complete derivation of the optimal broadcasting maps in the universal case and 
in the phase covariant case respectively. In Sect. |V|we study the entanglement properties of the states of two copies 
at the output of the universal and the phase covariant broadcasting maps. Finally, in Sect. IVII we summarize and 
comment the main results of this paper. At the end of the paper we report the details of the calculations needed to 
derive the results presented in three appendices. 

II. SYMMETRIC QUBITS BROADCASTING 

In this Section we introduce in a unified theoretical framework some preliminary concepts that will be employed to 
describe covariant symmetric qubit broadcasting maps. These concepts will be then specified to the universal and the 
phase covariant cases in the subsequent sections. A main tool we will extensively use in deriving the optimal maps 
is the formalism of the Choi-Jamiolkowski isomorphism j2jj between completely positive (CP) maps £ from states 
on the Hilbert space ^ to states on the Hilbert space J%f , and positive bipartite operators R on J(f ® Jff '. Such an 
isomorphism can be specified as follows 

R = £®l(\n)(n\) «—►£(/>) =Ttjr[l®P T R], (1) 

where Q is the non normalized maximally entangled state ^ m \ip m ) ® \ipm) m ^ ® ^ i ^ gives the identity trans- 
formation and X T denotes the transposition of the operator X on the same basis \ip m ) used in the definition of 

n. 

In terms of R, the trace-preservation condition for the map £ reads 

It jf[R] = lje . (2) 

Suppose that the map £ is covariant under the action of a group G. In this case the covariance property is reflected 
to the form of the operator R by the following correspondence 

£(U gP Ul) ee V g £{p)Vl ^ [V g <8> U* g ,R] = 0. (3) 

In the above expression U g and V g are the unitary representations of G 9 g on the input and output spaces respectively, 
while X* denotes complex conjugation on the fixed basis \ip m }- In this framework it is also possible to study group- 
invariance properties of the map £ in terms of the operator R. In this case we have the following equivalences 

£(U gP Ul) ee £{p) [1 ® U* g ,R) = 0, (4) 

and 

£{p) ee V g £{ P )V^ g [V g ® 1, R] = 0. (5) 

The above expressions refer to invariance properties on the input and output spaces respectively. 

In the following we will consider maps B from states of N qubits to states of M qubits, namely CP maps from states 
on ffl — (C 2 ) (8JV to states on J(f = (C 2 )® M . We will consider in particular symmetric broadcasting maps, namely 
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transformations where all receivers get the same reduced state. The figures of merit which are commonly used are 
invariant under permutations of the output copies, and this allows to assume that the output state of a broadcasting 
map is permutation invariant without loss of generality. Moreover, since the input consists in N copies of the same 
state, there is no loss of generality in requiring that the map is also invariant under permutations of the input copies. 
These two properties, according to Eqs. l@J|5j), can be recast as follows 

\n?®H?,R] =0, Va,r, (6) 

where 11^ and 11^ are representations of the output and input copies permutations, respectively. Notice that per- 
mutations representations are all real, hence II* = 11^. 

A useful tool to deal with unitary group representations U g of a group G on a Hilbert space Jif is the Wedderburn 
decomposition of J4? 

JT ~ ® , (7) 

where the index \x labels the equivalence classes of irreducible representations which appear in the decomposition of 
U g . The spaces Jf^ support the irreducible representations and C dfl are the multiplicity spaces, with dimension 
equal to the degeneracy of the /x-th irreducible representation. Correspondingly the representation U g decomposes as 



^ = 0^®^. (8) 

where 1^ is shorthand for Ij-.^ . By Schur's Lemma, every operator X commuting with the representation U g in 
turn decomposes as 

In the case of permutation invariance, the so-called Schur-Weyl duality [2l| holds, namely the spaces C d f for 
permutations of M qubits coincide with the spaces for the representation Uf M of SU(2), where U g is the defining 
representation. In other words, a permutation invariant operator Y can act non trivially only on the spaces J^Lj 
namely it can be decomposed as 

F = 0F M ®l dM . (10) 

The Clebsch-Gordan series for the defining representation of §U(2) is well-known in the literature (see for example 
pH |22| ) , and its Wedderburn decomposition is given by 

M 1 2 

3V ~ ^5 ® C d i , (11) 

0= do 

where Jfflj = C 2j+1 , jo equals for M even, 1/2 for M odd, and 

dj = + 1 ( M V (12) 

3 M/2+j + l\M/2-jJ V > 

In the case of the broadcasting maps the Hilbert space Jif % ffl on which the operator R acts, supports the 
two permutation representations corresponding to the output and input qubits permutations. Therefore it can be 
decomposed as 

/M/2 \ I N/2 \ 

0^®C d M® 0^®^' . (13) 

\j=jo J V=/o / 

By rearranging the factors in the above expression, we can recast the decomposition in a more suitable way, namely 

M/2 N/2 

jtr®.jf ~00(^(g).^?)(g)(c tij ®c d! ) . (w) 

j=jo l=lo 
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The operator R, in order to satisfy the permutation invariance property JBJ, according to Eq. (|10(l . can be written in 
the following form 

M/2 N/2 

R = ©0^® (H , (15) 

j=jo l=lo 

where the operators Rji act on J#j Cg> Moreover, in order to fulfill the requirements of trace preservation and 
complete positivity, the operators Rji must satisfy the constraints 

RflZO, Tr^ j7 ] = ^±i (16) 
a j 

where Tij denotes the partial trace performed over the space Jffj in the j-th term of the decomposition in Eq. 114J1 . 
and I2/+1 is shorthand for 1,^. We have now all the tools to study symmetric qubits broadcasting devices. In this 
work we are interested in the case of covariant broadcasting maps, which in general have to fulfill also the following 
covariance condition under the representations V? and V® M of a group G (see Eq. @) 



v® m ®v;® n ,r 



0. (17) 



The above condition gives a further constraint on the form of the operators Rji in Eq. (|15f) . Actually, the group V g is 
in general just a subgroup of the defining representation U g of SU(2), and therefore the representation V® M ® V*® 
acts non trivially only on the subspaces ® which are the ones supporting the operators Rji. In the next 

sections we will consider two interesting cases, namely V g = U g and V g = = e l ^ u " , corresponding to universal and 
phase covariant broadcasting respectively, and we will see how the form of the operators Rji depends on the particular 
choice of the considered covariance group. 

In addition to the Wedderburn decomposition and the related Schur-Weyl duality reviewed above, another useful 
tool we will extensively use in the following is a convenient decomposition of an iV-partite state of the form p® N , 
representing N qubits all prepared in the same generic state p. In Appendix^ we report the complete derivation of 
the following identity, which was originally presented in [T?) , 

N/2 j , 

^ = (r + r_f /2 E (Jr) \j™)(jrn\®t dr (18) 



3=]o m=—} 



For the sake of simplicity, in the above expression we considered density operators of the form p = (U + ro~ z )/2, 
r± = (1 ± r)/2, namely qubit states whose Bloch vector (of length r) is aligned along the z axis, and consequently 
the states \jm) are eigenstates of the operator J z in the j representation, namely J^P = Yj* m =-j m\jm)(jm\. Notice 
actually that the total angular momentum component J z of N qubits is clearly permutation invariant and therefore 
it can be written as 

N/2 

J z =0jW®l d .. (19) 

j=jo 

We want to point out that the decomposition (|18fl holds for any direction of the Bloch vector, provided that the 
eigenstates of J z in Eq. I|18l) are replaced by the eigenvectors of the angular momentum component along the 
direction of the Bloch vector in the single qubit state p. 

We will prove in Appendix C that the single-site output copy p' of a covariant broadcasting map commutes with the 
input density operator p. In order to quantify the performance of the broadcasting map B and to judge the quality 
of the single-site output density operator p' = Trjv/_i [B(p® N )] we will then evaluate the length r' of its Bloch vector, 
namely 

Tr[a zP '} = r' . (20) 

Notice, moreover, that the length of a Bloch vector r is simply related to the purity Tr[p 2 ] of the density operator 
p as Tr[p 2 ] — (1 + r 2 )/2. Therefore, maximizing the output Bloch vector length r' is equivalent to maximizing the 
output single-site purity. Notice also that so far we cannot exclude that the input and output Bloch vectors r and r' 
are antiparallel, and this just implies that r' can range from —1 to 1. 
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We will now show how to evaluate r' according to Eq. i|20|) . which is the main quantity that we will consider in 
the next sections in the particular cases of universal and phase covariant broadcasting. The trace in Eq. I|2(JI) can be 
computed by considering that the global output state S = B(p® N ) of the M copies is by construction invariant under 
permutations, hence 



r> = Tr[a z p>] = Tr [(a z ® 1® M ~ 1 ) S] = ± Tr 



]TlI CT (tr,®!®^- 1 )^ 



(21) 



The last term on the r.h.s. of Eq. Ij21|l contains a sum over the Ml possible permutations of the M output qubits. 
Notice that 



1 1 M 9 

-iVn.fe®!®^- 1 ))^ = — Y of = — J z , (22) 

where the operator a z acts as a z on the z-th qubit and identically on the remaining qubits. Now, by exploiting the 
permutation invariance of E, we can write 

M/2 

E = E i ® t dj , (23) 

3=3o 

and clearly 

2 M/2 

^-gE^^^l' ( 24 ) 

3=30 

The explicit expression of r' in the universal case will be derived in Sect. II I II In the phase covariant case, we will see 
that it is more convenient to take p diagonal on the a x eigenstates. The above formula in this case is just substituted 
by 

M/2 

r ' = ME^ Tr [ J ^-] (25) 

3=30 

and will be explicitly calculated in Sect. IIVI 

We want to stress that maximization of the figure of merit r' allows to optimize the fidelity criterion as well. In 
fact, for r' < r the two criteria coincide, whereas for r' ^ r one can always achieve unit fidelity by suitably mixing 
the output state with optimal r' and the maximally mixed one. On the other hand, direct maximization of fidelity is 
not analytically feasible, since fidelity is a concave function over the convex set of covariant maps, whence it is not 
maximized by extremal maps. 

Finally, we want to mention that in the next sections we explicitly maximize the scaling factor for N inputs and 
M outputs p N ' M (r) = r 1 It, which can be referred to as shrinking factor or stretching factor, depending whether 
it is smaller or greater than 1, respectively. It is obvious that this maximization is equivalent to maximizing r'. 
Superbroadcasting corresponds to the cases where p(r) > 1. 



III. UNIVERSAL CASE 



In this Section we will give the explicit derivation of the optimal universal broadcasting maps. Starting from the 
general broadcasting map described in the previous section we have to impose in this case the additional constraint 

[U? M ®U;® N ,R] =0, (26) 

where U g is the defining representation of the group §U(2). For the defining representation U g the following property 
holds 



U g = O-yUgCTy 



(27) 
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By exploiting such a property, the commutation relation (|26|l can be written more conveniently as follows 

[U^ M+N \S] = 0, (28) 

where S = (1® M <g> a® N )R(t® M ® af N ) = (t® M ® e l7rJ y) R (t mi ® e""^). The complete positivity and trace- 
preservation constraints in terms of the operator 5* are then equivalent to 

S^O, Tc J r[S\ = l je . (29) 
Upon defining Sji — ( i-2j+i ® e l7rJ y ' J i?^; (l 2 j+i ® e _l7rJ « ' ) , the constraints for complete positivity and trace preser- 



vation are now given by the following conditions on the operators Sji 



S S i>0, Tr j [S jl } = ^. (30) 



By exploiting the fact that the Clebsch-Gordan series for Jfj ® J%i is just ©jt-i^ji ^j, we can write 

M/2 N/2 j+l 

<8 = J^'' ! ® C d ' <g> C d < . (31) 
i=jo '='o J=|i-J| 

Notice that this is not the Wedderburn decomposition, since not all the subspaces J&j' ~ C 2,7+1 support inequivalent 
representations. However, the Wedderburn decomposition can be recovered by a suitable rearrangement that takes into 
account the repetitions of the same representation J. Using the decomposition (|31|) we can formulate the constraint 
(|28|l in terms of the operators Sji as follows 

j+l 

S H = 4l P f,i> ( 32 ) 
J=\j-l\ 

where, by complete positivity, the coefficients \ are real and positive, and P^ is the projection of the space 3#j®J%i 
onto the J representation, satisfying 

r Ti 2J+1 r M 2J+1 

T ^ [Pf,l] = 2/+T laI+1 ' Tn [jR " ] = 2j + T l5tf+1 ' (33) 
The set of projectors Pj 7 ; is clearly orthogonal. The trace-preservation constraint (J5Jl can now be written as 

N/2 M/2 j+l o T , -i 

0E E ^ 1 ^ r d i (l 2m ®l d! ) = l^ ! (34) 

i=io J=J0 J=|i-i| 



which is equivalent to the conditions 

M/2 j+l 

21 + 1 



E E 1 (35) 

j=jo J=\j-l\ 

Along with the complete positivity constraint sj \ ^ 0, Eq. (|35() defines a convex polyhedron whose extremal points 
are classified by functions j = ji and J — J\ 

j 2? + 1 1 

^ 2T^r^ (36) 

The classification of symmetric universally covariant maps is then completely determined in terms of the vectors ji and 
Ji, whose elements can range from jo to M/2 and from \ ji — l\ to ji + 1, respectively. Extremal maps then correspond 
to the following form for the operators S 

N/2 2/ + 1 1 

s-Qur+i^* 1 **® 1 *- (37) 
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The optimization of the figure of merit r' or, equivalently, of the scaling factor p(r) can be obtained by explicit 
calculation from Eq. (|24H . The output state £ of the broadcasting map applied to an input state p® N can be 
represented as 



S = Tr^r[l® M <g> (a y p T a y )® N S] = Tr.*.[l® M <g p® N S] 



(38) 



where p denotes the orthogonal complement of p, which just corresponds to the change r — > —r (or, equivalently, 
r± — > Tzf). Using the decomposition in Eq. (|18|l for p® N and the form l|37[l for the operator S, we can express Eq. (|38[) 
as follows 



N/2 



Z=Z Jl n=-l v 

We can now use Eq. (|24|l to evaluate the scaling factor, namely 



Tr, 



(t 2jl+1 ®\ln)(ln\) Pp a 



(39) 



2 , , .%•/•> ^ — •> 21 + 1 



^ W Mr l+ 1 f-f 2 J; + 1 ' ^ 



/— /o n— — l 

N,M 



Tr 



(40) 



In Appendix FBI we report the explicit calculation of p N - M (r) and we show that it can be written in the following form 



N/2 



where 



J r (J r +l)-j(j + l)-/(f+l) 

2Z(Z + 1) 



(41) 



(42) 



Since r_ r+, the sum _, n(r_ /r+) n in Eq. 141(1 is always negative. Therefore, the function p N ' M (r) is 
maximized by the choice of J; and ji minimizing (3, which clearly implies Ji — \ji — l\. The form of the coefficient 
P(Jl,jl,l) for ji < I is given by 



ji,ji,l) = -j , 



whereas for ji /? Z we have 



31 + 1 
l + l 



(43) 



(44) 



In both cases (3 is minimized by choosing the maximum value of ji, and therefore the maximum scaling factor is 
achieved by ji — M/2. For M > N the optimal value of the figure of merit is then univocally determined by the value 
of the function (3 



f3(M/2-l,M/2,l) 



M + 2 
2(1 + 1) 



while the optimal scaling factor is given by 

» N > M ( r \ = M + 2 
Mr 



p"> M (r) 



N/2 I 

di 



l=lo 



n——l 



(45) 



(46) 



The corresponding output state takes the form 

N/2 



l=l 



M — 21 + 1 



n=-l m=-M/2 



M , . . , _ 

— l,m + n ) — 



M \ I M 
2 ' / \ 2 ' 



(47) 
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p (r) 




FIG. 1: Optimal scaling factor p N ' N+1 (r) = r' / r versus r for the universal broadcasting, for M = N + 1 and N ranging 
from 10 to 100 in steps of 10. Notice that there is a wide range of values of r such that p N ' N+1 (r) > 1, corresponding to 
superbroadcasting. 



where (4r-m, ln\^- — l,m + n) denote the Clebsch-Gordan coefficients. 

As mentioned in the previous section, in Appendix Icl we prove that the single-site reduced output state Tr&f_i[E] 
commutes with a z , hence p N,M (r) is definitely a scaling factor. Two interesting cases we will consider in the following 
are the ones with M — N + 1 and M — oo, for which the scaling factor takes the explicit forms 



*r, o N / 2 A I / 

P [ ) (N +l)r lr+r - j 2^1 + 1 2^ n [ r+ 

x t— to n— — I x 

1 N/2 A 1 / \» 

l=l n=-l V + ' 



(48) 



The function P N ' N+I (r) is plotted in Fig. [T]for N ranging from 10 to 100 in steps of 10. We can see that for a suitable 
range of values of r the scaling factor is larger than one. This corresponds to a broadcasting process with an increased 
single-site purity at the output with respect to the input. This phenomenon occurs for N ^ 4. In this case p(r) is 
actually a stretching factor, and we call such a phenomenon superbroadcasting. The maximum value of r such that 
it is possible to achieve superbroadcasting will be referred to as r* (N, M) and it is a solution of the equation 

p N > M {r*) = 1. (49) 

It is clear that the optimal scaling factor for fixed N is a non increasing function of M. Actually, by contradiction, 
suppose that the map with M + K output copies has a higher purity than the optimal map with M copies. Then one 
could trace over K copies from the former map, and he would obtain a map with M output copies with purity higher 
than the optimal, which is obviously absurd. This implies that in general r*(N,M) < r*(iV, M + K), and for large 
values of K superbroadcasting may not be possible anymore. The maximum M such that superbroadcasting can be 
achieved for N input copies will be referred to as M*(N). It turns out that, apart from the values N = 4, 5 for which 
we have M*(4) = 7 and M*(5) = 21, for TV ^ 6 one has M*(N) = oo, namely superbroadcasting is possible for any 
number of output copies. In Fig. [5] we report the values of 1 — r*(N,N + 1) and 1 — r*(N, M*(iV)) for 4 ^ TV ^ 100. 
By a numerical analysis we have evaluated the power laws for the two curves, which turn out to be in good agreement 
with 2/7V 2 and 1/N, respectively. 

We want to point out that for input pure states (r = 1) only the term with I = N/2 in the expression (|37fl is 
significant. The optimal map then corresponds to the optimal universal cloning for pure states derived in |l4j . 



IV. PHASE COVARIANT CASE 



In this section we study the case of symmetric phase covariant broadcasting, where we restrict our attention to input 
states lying on an equator of the Bloch sphere, say the sy-plane. The equatorial qubit density operator in this case has 
the explicit form p = (1 + r cos 4>a x + r sin 4>a v )/2. The starting point, as in the case of universal broadcasting, is the 
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FIG. 2: The logarithmic plot reports the behaviour in the universal case of 1 — r*(N,N + 1) and 1 — r*(N, A'I*(N)) for 
4 ^ N ^ 100. The upper line, corresponding to M = M t (N), has a power law 1/N. The lower line, corresponding to 
M = N + 1, has a power law 2/N 2 . 

requirement of permutation invariance © for input and output copies, that leads to the form l|15fl of the operator R. 
Moreover, in this case we demand covariance under the action of the group of rotations along the z-axis = e^ CT */ 2 . 
We proceed analogously to the case of universal broadcasting. By imposing the covariance condition for the map we 
require invariance of the R operator, namely [V® M ® V^® ,R] = 0. By exploiting the Wedderburn decomposition 

© for the operator V® M , namely 

M/2 

Vf M = 0^" ®t dj , (50) 

3= jo 

the phase covariance requirement corresponds to the following additional condition for the operators Rji 

fe.e^^e-^Vo, Vj,I, (51) 
where is defined according to Eq. (|19|l . A convenient form for the operators Rji satisfying Eq. (|51Jl is the following 

l i j-i 

%=EE E ^ n , tk \j,n + k)( 3 ,n> + k\®\l,n)(l,ni\ (52) 

n— — l n' — — l k—l—j 

when j ^ Z, and 

R oi= E E E ril tm , tk \j,m)(j,rn'\®\l,m + k)(l,m' + k\ (53) 

m——j m' — ~j k—j — l 

when j < I. Notice that there arc two more running indices with respect to the universal case (|32ll . The index 
n' in Eq. 1)521) simply allows for off-diagonal contributions in the operator Rji, while we will see that the index k, 
which labels equivalence classes, is related to the direction of the reduced output state Bloch vector. In particular we 
will show that, in order to get an equatorial output, the operators Rji have to be symmetric in fc, in the sense that 

r n n' k — r nn' -k- Notice also that k takes integer values when M — N is even and half integer values when M — N 
is odd. 

The trace-preservation condition l|lb[l now reads 

M/2 

E E <n,fcdi = l. VZ,«, (54) 

j=jo k=~\l-j\ 

and, analogously to the universal case, the fact that the operators Rji are diagonal with respect to the indices j's and 
/c's implies that the extremal points are classified by functions 




j=ji, k = ki, 



(55) 
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and satisfy 



n,n : ki 



(56) 



We will now compute the output density operator and the scaling factor for N — > M phase covariant broadcasting 
maps. Without loss of generality, let us now consider an input state p oriented along the x-axis, namely p = (Tl+ra x )/2. 
The density operator p® N can then be decomposed, analogously to Eq. JTSJ, as 



N/2 I 

©I 

I—Iq n— — l 



\r, n )(r,n\ <g> l dl 



(57) 



where \l x , n) is the eigenvector of Ji corresponding to the eigenvalue n. In the following, eigenvectors without explicit 
specification of the superscript axis, such as \jm), are intended to be along the z-axis, namely \j z , m). According to 
Eq. QJ, the density operator Son Jfs (C 2 )® M , describing the output state of the M copies, can be written as 



£ = Tr^ [(1® M ® (p T ) 8JV ) i2] = Tr^ [(1® M ® R] 



M/2 

j—jo l.k n,n' n" 



(Wi) 

\j, n + k) (j, n' + k\ <g> Id, i 



(58) 



where (Wfc) & = (fe, aj A; 31 , 6) are the entries of the Wigner rotation matrix in the k representation which rotates the 
^-components into the x-components — in the usual notation (the one that is found, for example, in |22() such entries 

are denoted as (/3 = -|). As discussed previously in Sect.[nJ the projection r' along the x axis (|25|1 of the Bloch 
vector of the single-site output state is a convex (linear) function on the convex set of phase covariant broadcasting 
maps, and therefore it achieves its maximum on extremal broadcasting maps. Let the functions j = ji and k = k\ 
denote an extremal map. Hence, starting from Eq. (|25|l and specializing Eq. I|B11|I . derived in Appendix B, to the 
extremal case j = ji and k = k\ , we can express the scaling factor in the following form 



JV.M 



N/2 
\N/2ST 



I 



;^; n——l 



Jul 



exp 



lof 



1 



n,n+l 



j(h) 



n+ki ,n+fei + l 



(59) 



where [X^} n ^ m denotes the matrix element of the operator evaluated with respect to the eigenstates of i. e. 

\jm). The final form in Eq. I|59(l is now suitable to be optimized. First of all, since the matrix elements j| are 

non-negative, the maximum purity is reached by maximizing the off-diagonal elements of Rji, namely for rank-one 
Rji with all the matrix elements equal to 1/dj, (see Eq. |(SSJ). We now want to identify the values of ji and fc; 

corresponding to the optimal scaling factor of the map. The matrix elements of Jx 1 ^ take the explicit form 



jUi) 



1 



t , , = - VllUl + 1) - (« + ki){n + h + l) 

n-\-ki ,n-\-Ki +1 ^ 



(60) 



Since the above matrix elements are maximized in the central block of the matrix, the optimal map is achieved by 
choosing ki as close as possible to zero, for all the values of I. When M — N is even, the optimal choice corresponds 
to ki = for all /. When M — N is odd there are two equivalent possible choices for ki, namely ki = ±1/2, for each 
value of /. Moreover, we set ji — M/2 for all I, namely as large as possible. For M — N even, the global output S 
and the scaling factor are given by 



exp ( 4° log 



N/2 I 

1=Iq n 7 n'——l 

A N/2 l-l 

l—lo n— — l 



exp 



1 + r 
1 - r 



A l hot 



M \ I M . 
2 A 2 



1 + r 
1 -r 



(61) 



n,n+l 



j(M/2) 



n,n+l 



For M — N odd we have many more solutions, corresponding to all the possible combinations of fc; = ±1/2 for all 
values of /. As will be clear in the following discussion, we will examine the two cases of fe; = 1/2 and ki = —1/2 for 
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all I. In the former case we can write 



N/2 i 
I—Iq n.7i / — — l 

a N/2 r 

I—Iq n— —I 



l+r 
l-r 



J« lot 



M 1\ M . 1 



(62) 



l + r 
l-r 



j(M/2 



while for fc; = —1/2 we have 



N/2 I 

Z=/ n n.n' — — l 
. N/2 l-l 



P'o 



Notice that, since 



r(j) 



m,m+l 



Z=/o n— — l 



exp 



1 + rV 




M 


1-rJ. 




Y 


(4" log 


i + 




l-r/ 



1\ / M 



n,n+l 



j(M/2) 



n+l/2,ri+3/2 



ri-l/2,n+l/2 



(63) 



T (J) 



-(m+1), — m 



and the same property holds for the matrix elements of any power 



of J Xl the scaling factors corresponding to the extremal maps with ki = 1/2 and ki = —1/2 are exactly the same. 
This means that the Bloch vector components in the xy-pl&ne are scaled in the same way by the two maps. 

We want to point out that an extremal map with ki generates output density operators with a non vanishing 
component of the Bloch vector along the z direction. Actually, for the input state (|57|) the output single-site density 
operator is given by 



where 



p' = ~(1 + r'a x + a k a z 



(64) 



(65) 



Optimal broadcasting maps, where the Bloch vector is just scaled along its input direction, can then be obtained for 
odd values of M — N by equally mixing the two maps considered above, corresponding to ki = 1/2 and fc; = —1/2. 
As mentioned earlier, since the two maps give the same scaling factor, their mixture does not compromise optimality. 
Notice that the optimal broadcasting maps we have derived in this way are independent of the input state. In the 
limit of pure input states, the above maps coincide with the optimal phase covariant cloning for pure equatorial states 
presented in Ref. p^ |. 

We will now discuss more quantitatively the results derived above. The optimal scaling factors, reported in Eqs. 161|) 
and (|62|) . contain only known terms and can be studied numerically. It turns out that phase covariant superbroad- 
casting is possible even for N = 3, with M*(3) = 12. Moreover, it is possible to superbroadcast an infinite number 
of output copies starting from N = 4 (M*(iV) = oo for N > 4). As for the universal case, we can easily compute 
the function p{r) for M = N + 1 and M = oo, which is monotone decreasing in M. In Fig. |3|we report the plots of 
p Ar,Ar+1 (r) for values of N such that 4 < N < 100 in steps of 8, and for M = N + 1. In Fig.0]we report the plots of 
the values of 1 — r*(N,N + 1) and 1 — r*(N, M*(N)), as defined in the universal case. The upper line refers to the 
case N — > oo and shows a behaviour like 1/2N. The lower line is for N — > N + 1 and scales like 2/(3N 2 ). 

As before, in the limit of pure input states (r = 1), the optimal phase covariant superbroadcasting map coincides 
with the optimal phase covariant doner for qubits of Ref. |24j . 



V. BIPARTITE ENTANGLEMENT IN THE GLOBAL OUTPUT STATE 



In this section we analyze the entanglement properties of the output state of the optimal N — > M broadcasting 
maps. Notice that, since broadcasting maps are always optimized with ji = M/2, the output state is supported on 
(which has multiplicity d/vf/2 = 1); namely the completely symmetric subspace of (C 2 )® M . Therefore, also the 
reduced state of two qubits = TrM-2p] is symmetric. We will analyze in particular bipartite entanglement in 
the output state, which is conveniently described in terms of the concurrence [26j 



C{p {2) ) = max{0, Ax - A 2 - A 3 - A 4 } , 



(66) 
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FIG. 3: Optimal scaling factor p N ' N+1 (r) = r'/r versus r for the phase covariant broadcasting, for M = N + 1 and N ranging 
from 4 to 100 in steps of 8. 
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FIG. 4: The logarithmic plot reports the behaviour of 1 — r*(7V, N + l) and 1 — r,(N, M* (N)) in the phase covariant broadcasting, 
for 3 ^ TV ^ 100. The upper line, corresponding to M = M,(N), has a power law 1/2N. The lower line, corresponding to 
M = N + 1, has a power law 2/3N 2 . 



where A, are the decreasingly-ordered eigenvalues of the operator $ = y \J p^p^y/ p( 2 \ and p^ — er® 2 /?( 2 )*er® 2 . 

We will first consider the universal case, where the output state S is diagonal on the J z basis. As shown in 
AppendixEl the state p^ commutes with and therefore it can be written as a linear combination of independent 
powers of namely 

p& = a l+/3ji 1 )+ 7 (ji 1 )) 2 . (67) 

In the above expression the positivity and unit trace constraints are given by 

a + 7^|/3|, cO0, 3a + 2 7 = l. (68) 

The eigenvalue A4 of \& is always 0, corresponding to the null component of p^ on the singlet. By the unit trace 
condition we can express 7 as a function of a and (3 as 7 — (1 — 3a) /2, and the positivity condition in terms of the 
two independent parameters (/?, a) is just 

a < 1-2101. (69) 

The above inequality defines a triangle with basis [—1/2, 1/2] and height [0, 1], as shown in Fig. [S] (left). A state in 
,3^1 is then completely determined by the couple (/3, a). Notice that the only pure states of the form (|67|l are |1, 1), 
|1, 0) and |1, —1), which correspond to the vertices (1/2, 0), (0, 1) and (—1/2, 0) respectively of the triangle in Fig. 
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We will now express the concurrence in terms of (/?, a). Since p^ is real it follows that p^ = a® 2 p^ a® 2 , and 
therefore we can write 

=e i * J PpWe- i * J P. (70) 

It is easy to verify from Eq. (|67(l that p^> corresponds to the couple (a, —0). Moreover, since p^ and p^ commute, 
the operator 'J can be simply written as \J p( 2 )p( 2 ). By exploiting some algebra, and taking into account the identities 

\Jzj = and (j* 1 ^ = (j^J ' we S et tne following expression 

p^pW =a 2 t + (2a 1 + 1 2 -f3 2 )(jpy . (71) 
^From the above expression, by using the unit trace constraint, we can compute the eigenvalues of ^ 



l ^-^-^- ^j (72) 

Notice that the first eigenvalue is doubly degenerate. The concurrence can then be written as follows 

1-4/3 2 

> - y/l - 2a + a 2 - 4/? 2 , a> 

The above equation defines a parabola inside the triangle l|69|) of states. Such a parabola separates the region of 
separable states from that of entangled states, shown in light and dark gray in Fig.Elrespectively. In order to analyze 
the amount of bipartite entanglement in the broadcast states, we have then to evaluate the couple (/3, a) for the 
reduced state of two output copies and then determine in which region of the triangle it lies. Using Eq. (jUTo| derived 
in Appendix C, we can numerically evaluate (/3, a) for the universal double-site reduced output density operator 
pi 2 ) — TrM-2p]- In Fig. 0we report the parametric plot for the case of 4 input and 5 output copies. As we can see, 
the black line moves towards positive (3 as the Bloch vector length r of the input state goes from to 1. It is possible 
to see in the magnified plot on the right that, as r gets close to 1, i. e. the input state gets pure, the output exhibits 
bipartite entanglement, since it crosses the parabola. In the limit of pure input states, these results agree with the 
ones derived in Ref. |27j . 




FIG. 5: The figure represents the set of bipartite symmetric states which are diagonal on the jjft basis. All such states are 
parametrized by the couple (/3, a), as given in Eqs. I67H and The light gray region contains separable states. The dark 

gray region contains the entangled states. The black line is the parametric plot of the double-site reduced output of the 4 — »■ 5 
universal broadcasting, for input Bloch vector length r ranging from to 1. In the magnified window on the right, it is possible 
to notice that, for nearly pure input, the black line crosses the parabola, namely, the output exhibits bipartite entanglement. 



In the phase covariant case it is not possible to carry on the same analysis, since, as we notice in Appendix El the 
global output state does not commute with J x . However, using the partial traces in Eqs. IjClOl) . I|C11|) . and iC12f) . it 
is still possible to evaluate the concurrence numerically. In Fig. we report the plots of the entanglement E defined 
in Ref. 12a as follows 



l + Vi - c 2 " l + yi - W l-VT - c* i-VT~c^ 

E= = log log . (74) 
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FIG. 6: Entanglement E for the double-site reduced output state, in the cases of universal (l.h.s) and phase covariant (r.h.s) 
broadcasting from N = 2, 4, 6, 8, 10 to M = N + 1 copies, as a function of the input Bloch vector length r. 



for N = 2, 4, 6, 8, 10 and M = N + 1 as a function of the input Bloch vector length r, both in the universal and phase 
covariant case. Notice that, contrarily to what happens in the universal case, in the phase covariant case bipartite 
entanglement vanishes in the limit of pure input states. The absolute value of C goes to zero for increasing number 
N of input copies. 



VI. CONCLUSIONS AND FURTHER DEVELOPMENTS 



In this paper we studied symmetric broadcasting maps, where N input qubits initially prepared in the same mixed 
state are transformed into an output state of M qubits, all described by the same density operator. We considered 
covariant maps and we investigated in particular the universally covariant case and the phase covariant case. We have 
shown that for sufficiently mixed initial states and for N ^ 3 it is also possible to partially purify the single qubit 
output density operator in the broadcasting operation. Such a new process was named superbroadcasting. 

The new superbroadcasting channels open numerous interesting theoretical problems. The first problem is to extend 
the map to any dimension d > 2 and to different covariance groups. Indeed, for special cases it is easy to see that 
increasing the dimensionality and/or reducing the set of states to be broadcast makes superbroadcasting possible 
with smaller N, and even with N = 2 input states. As a matter of fact, this is the case of universally covariant 
superbroadcasting from N — 4 to M — 6, which can also be regarded as superbroadcasting for d = 4 for special states 
of the form p ® p, and for the covariance group SU(2) x §U(2). The case of dimension d — 4 is most interesting, 
since it can be exploited to improve entanglement for bipartite states of qubits. Also the infinite dimensional case 
(the so-called "continuous variables" ) turns out to be interesting, and easily feasible experimentally . It should be 
emphasized that for dimension d > 3 there are many ways of increasing purity, and certainly the most interesting case 
is the purification along the mixing direction of a noisy channel (notice that most channels do not correlate, whence 
the produced state is the tensor product of identical mixed states) . 

Another major problem is the analysis of the detrimental correlations between two outputs, e. g. to establish 
whether they are quantum or classical. These correlations are exotic, in the sense that instead of increasing the local 
mixing as usual, they reduce it. Such mechanism is new, and deserves a more thorough analysis. An interesting 
issue, for example, is that they cannot be erased leaving the local state unchanged (the de-correlating map — which 
sends a state to the tensor product of its partial traces — is non linear) , and this raises the problem of the optimal de- 
correlating channel, which optimizes the fidelity between the input and the output local state. Such optimal channel 
can be derived using the same technique for optimal covariant maps used in the present paper. 

Finally, distributing quantum information — and in particular superbroadcasting — raises the new problem of the 
trade-off between broadcasting and cryptographic security. Indeed, on one side, the presence of many identical uses 
seems to open more possibilities of eavesdropping, however the detrimental correlations may drastically reduce such 
possibility, and the opportunity of detecting the eavesdropping on the joint output state may be exploited to increase 
the security. 
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APPENDIX A: DECOMPOSITION OF p® N 

The global density operator p® N is clearly invariant under permutations of the N qubits, and, according to the 
Schur-Weyl duality, it can be represented by the Wedderburn decomposition I|1U|) 

N/2 

P W = ©Pi®l^ (Al) 

3= 3a 

where pj is a state on jtfj. In order to evaluate pj it is sufficient to evaluate the matrix elements 
(jm\ £g> (a\p® N \jm') <g> |a), where \jm) are eigenstates of J z in the j representation, and \a) is an arbitrary state 
in C dj . For the sake of simplicity we can suppose that the state p has the Bloch form i(l + ra z ). The problem is 
now to choose \a) in a suitable way. It turns out that a clever choice is given by 

\jm) <g> |a) = \jm) 2j ® l* - )®^, (A2) 

where the subscript 2j means that \jm)2j is a vector in the symmetric subspace of the first 2j qubits, while \^~) 
is a singlet, supporting an invariant representation on a couple of qubit spaces. Notice that in Eq. (|A2|) the tensor 
product on the l.h.s. refers to the "abstract" subspace Jtfj (B> C dj in the Wedderburn decomposition, whereas the 
one on the r.h.s. refers to the decomposition (C 2 )® 2 - 3 ® (C 2 )®^ -2 - 7 ') grouping separately the first 2j qubits and the 
remaining N—2j. Moreover, since the chosen density operator p® N commutes with the total J z , given by the following 
expression 

N/2 N/2 j 

J* = 4 3) ® H = ]T m\jm){jm\ ® l dj , (A3) 
3=3o j=jo m=—j 

then pj commutes with , which implies that pj is diagonal on the eigenstates | jm) of J z ^ . Therefore we can write 



( Pj )mm = viMp^ljmhj ((*- |p® 2 | *-))^" . (A4) 

Since |jto)2j is symmetric, it is a linear combination of factorized vectors with (j + m) qubits in the 1 1/2, 1/2) state 
and (j — m) in the state 1 1/2, —1/2). As a consequence we can also write 

p® 23 \jm) 2 j = r 3 + m r J S m \jm) 2 j , (A5) 

where r± = ^rp. By analogous arguments it follows that 

p® 2 \^-) =r + r_|* _ ). (A6) 

The matrix element ( P j) mm has the following expression 

(\ m 
. ( A7 ) 

and the decomposition of is finally given by 

N/2 j 



3=3o m=—j 



p ®N = {r+r £ — b'WO'H®^-- (as) 



Notice that this expression exhibits a singularity for r = 1 due to the rearrangement of terms (|A7|1 . However, a finite 
limit for r — > 1 exists, as it can be seen from the equivalent expression 

N/2 j 

P® N = (r+r_ ) ^ £ < +m rr m | jm) (jm| ® l d . , (A9) 

j=3o m=-j 

which exhibits no singularities. 
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APPENDIX B: FORMULAE FOR THE SCALING FACTORS 

In this Appendix we will derive the explicit form of the scaling factor for the universal and phase covariant cases. 



1. Universal case 



In order to calculate p N ' M (r) we start rewriting Eq. 14U|) as 



Tr 



(j® ® \ln)(ln\ 



P J 



Tr 



ln)(ln\Tr 3 [(jW®l 2 ,+i) P/, 



Let us define 



Xi l) = Tr, 



(j^ 5 ® lai+i) P/ 



P J 



where i = —1, 0, 1, and Xq = X z . Since 

(up*uw)p?,(up*uwy 

the set xf^ transforms according to 

UiOx^uf = Tr, [(t#>,l%W>t P/ ; 



= Tr, 



L \fc=-l 



(Bl) 
(B2) 

(B3) 



(B4) 



fc=-i 



and we conclude that is an irreducible tensor set. It can then be proved by the Wigner-Eckart theorem that 

xf 1 = ajf \ and in particular xjP = aji l \ From the last relation and from the identity 

\ (jf jW + J«jf ) + (J«) 2 = £ a,J«J« = (JW) 2 Z(i + 1) W, (B5) 

»=-l 

where a_i = ai = 1/2 and ao = a z = 1, and = we have 



*/(/ + 1)(2Z + 1) =a E a * Tr 

z=-l 
1 

= a 2 ; Tr 



j ( '])p j 



(B6) 



By using the well known identity 



i 

E a *( 



) Pf,l = 5 (^ (J)2 " J U)2 ® 12,+! - i 2j +l ® •/<'>*) P/ 



= P 



j J(J+l)-i(i + l)-J(I + l) 



we can write the explicit form of the coefficient a 

2 J + 1J(J + 1) - j(j + 1) - 1(1 + 1) _ 2J + 1 



2Z + 1 2/(2 + 1) 

By using the above expression we finally have 



2Z + 1 



Tr 



(B7) 



(B8) 



(B9) 
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Phase covariant case 



Substituting the global output state E given by Eq. (|55)l into Eq. (JStjJ, it is possible to compute the scaling factor 
p N < M (r) as 



i 



(BIO) 



n"=—l 



(W})n",n, 



where (Wk)ab = (k,a\k x ,b} are entries of the Wigner rotation matrix in the k representation which rotates the z- 
components to the ^-components (in the usual notation, the one that is found, for example, in |22| . such entries are 

denoted as (J3 = ^)). In Eq. (|B10|) the sum over n" gives the (n,rc.')-th matrix element of exp( log r + /r_). 
Therefore we can write 



P N ' M (r) = 77l( r + r -) Ar/2 E d J' d ' E 



Mr 

4 
Mr 



j,l,k n,n f — — l 



n— — / 



oxp ^ jj^ lof 
exp | lof 



1 + r 
1 -r 

1 + r 
1-r 



7-0) 



n,n+l 



n+/c,n+A;+l 



(Bll) 



where in the last equality we used the fact that J% = Yl m m \j X i m ) 0' X ' TO I nas non-null matrix elements only on the 
second-diagonals, and we multiplied the second line by a factor 2 considering in the sum only one of the two second- 
diagonals. Notice that all matrix elements in the previous equations are calculated with respect to the z-oriented 
basis \j z ,m) = \j,m) and \l z ,n) = \l,n). 



APPENDIX C: REDUCED OUTPUT STATES 



1. Single-site reduced output states 



In this appendix we want to calculate the following partial trace 

Tr M -i [\jm)(jm\®l dj ] , 



(CI) 



where the operator to be partially traced acts on Jtfj <g> C dj C (C 2 )® M , and \jm) are eigenstates of J* , as usual. In 
order to do that, we first decompose the vector \ jm) into its components onto j£j_i/ 2 ®J%i/2 using the Clebsch-Gordan 
coefficients 



\jm) = 



I j + m 
2j 



1 



1 



1 — i, Tn — — 
J 2' 2 



1 1 

2' 2 



U -m 
2j 



1 



1 



i i 

2' _ 2 



and then trace the operator \jm)(jm\ over Jfj_i/2- In this way we get 

1 m 



Trj-i/2 [\jm){jm\ 



2j 



(C2) 



(C3) 



We now recall a fact related to the already mentioned Schur-Weyl duality, by which multiplicity spaces C dj in the 
Wedderburn decomposition (jTT|l support irreducible representations of the permutation group {n£ 7 } of M qubits. 
Hence, for any operator O on Mj ® C dj one has 



Ml 



^nfonf = — tv, [o]®t dj . 



(C4) 



For convenience, let us write 



\jm){jm\ ® t d] = E 11 ^ (b' m >0'"*l ® I* )<*" 



(C5) 
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as we already did in Eq. (|A2|) . With this choice, we get: 

d, 



Ttm-i m){j, m\ ® t dj ] = Tr M 

d 



"2 
Ml 



(\j,m){j,m\® |MT><*~ 



n 



M 



-"J 

M\ 



1 



(M - 1)!(M - 2j)- + (M - l)!2j ( - + -a 



1 m 



(C6) 



The first term in the sum comes from excluding from trace one of the (M — 2j) qubits in singlet state. The second 
term in the sum comes from excluding from trace one of the 2j qubits in \jm) state, and from Eq. (|C3|I . Rearranging 
the above equation, we get the final expression 



Tr M _! [\j, m)(j, m\ <g> t dj \ = dj Q + j^v^j 



(C7) 



2. Properties of single-site output states 

Consider the global output state £ = B(p® N ). If it commutes with the total angular momentum component along 
the direction z, for example, then it is simple to prove that also p' — TrM-i[E] commutes with a z . In fact, from 
Eq. 1|C7(I . it is simple to see that 



oc a Z) (C8) 



Tr M -i [4 o) ® 1 

for all j, and consequently also Trj\/_i [J z ] oc a z . 

In the universal case, the global output S is diagonal on eigenstates of j| M ^ 2 \ hence p' = Tr^f-ip] commutes with 
cr z , according to previous arguments. In the phase covariant case, it is more difficult to prove on general grounds that 
[p' , <j x ] = 0, since in the phase covariant case [J x , S] ^ 0. The simplest thing we can do is to compute the partial trace 
of Eqs. (|61|l and (|62|l using again Clebsch-Gordan coefficients (IC2|) . First of all let us notice that, tracing over M — 1 
qubits, only the terms with |n— n'\ < 1 contribute. Among these, the terms with n — n' give one factor proportional 
to 1/2 and one factor proportional to a Zl whereas terms with \n — n'\ = 1 contribute with factors proportional to 

a x , since the matrix of coefficients exp ^ log is symmetric, see Eq. Ij64|l . Then, posing ctk — without loss of 
optimality, TrA/_i[S] commutes with a x . 



3. Double-site reduced output states 

We will show here how to compute the reduced output state of two copies, which is used in Section to compute 
the concurrence between two of the M clones. Clearly it does not matter which two clones we are considering, since 
the global output state is permutation invariant. Using Clebsch-Gordan coefficients, it is possible to decompose a 
vector in M'j into its components onto ® M\ 



\jm) =* r ; y ~ - 1.™ - i>im> + j^w-^ - 





+ m){j + 


m — 


1) 




2j(2j - 


1) 




(J 


- m)(j - 


m — 


1) 




2j(2j - 


1) 






and then to compute the partial trace of \jm)(jm'\ over For j = M/2 we have 

/ (M-2m)(7 

Tr M -2[|M/2,m)(M/2,m|] = 



(C9) 



/ (A-f-2m)(M-2m-2) n n 

' AM(M-l) U U 

\ n n (M+2m)(M+2m-2) 

\ U U 4M(M-1) 
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when |m — m'l = 1 



Tr A f_ 2 [\M/2,m)(M/2,m+l\ +h.c/ 
/ 







(M-2m~2)y/(M-2m)(M+2m+2) 
2V2M(M-1) 



(M-2m-2) A /(Af-2m)(M+2m+2) 
2V2M(M-1) 









\ 



V 

and |m — m'l = 2, 







(M+2m) A /(Af-2m)(M+2m+2) 
2 V / 2M(M-1) 



(M+2tn) A /(Af-2m)(Af+2m+2) 
2V2M(M-1) 





Tr A f- 2 [|M/2, m){M/2, m + 2| + h.c. 

/ 




n -y/(Af-2m)(Af+2m+4)(Af+2m+2)(Af-2m-2) \ 
U 4M(M-1) 





• v /(M-2m)(M+2m+4)(A/+2m+2)(A/-2m-2) „ 
\ 4M(A/-1) U 

For |m — m'| > 3, partial trace over M — 2 copies gives null contribution. 
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